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Abstract 

Chaotic properties of symmetrical two-dimensional stadium-like billiards with elliptical arcs are 
studied numerically and analytically. For the two-parameter truncated elliptical billiard the ex- 
istence and linear stability of several lowest-order periodic orbits are investigated in the full pa- 
rameter space. Poincare plots are computed and used for evaluation of the degree of chaoticity 
with the box-counting method. The limit of the fully chaotic behavior is identified with circular 
arcs. Above this limit, for flattened elliptical arcs, mixed dynamics with numerous stable elliptic 
islands is present, similarly as in the elliptical stadium billiards. Below this limit the full chaos 
extends over the whole region of elongated shapes and the existing orbits are either unstable or 
neutral. This is conspicuously different from the behavior in the elliptical stadium billiards, where 
the chaotic region is strictly bounded from both sides. To examine the mechanism of this difference, 
a generalization to a novel three-parameter family of boundary shapes is proposed and suggested 
for further evaluation. 
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I. INTRODUCTION 



Two-dimensional planar billiards are nonlinear systems with rich and interesting dynam- 
ical properties. A point particle, moving with constant velocity within a closed boundary 
and exhibiting specular reflections on the walls, can have regular, mixed or fully chaotic dy- 
namics, in strong dependence on details of the boundary shape. In physics, two-dimensional 
billiards offer good examples of coexistence of regular, mixed and chaotic dynamics in Hamil- 
tonian systems. This type of behavior, illustrated by the standard map and explained by 
means of the KAM-theorem, is present in many realistic phenomena, such as planetary sys- 
tems and various types of coupled oscillators [ij. Chaotic billiards were flrst introduced by 
Sinai who considered the defocusing effects of circular scatterers in the two-dimensional 
Lorentz gas. After the important discovery by Bunimovich 3|, |4| that also the focusing 
circular arcs can lead to a fully chaotic behavior, many investigations were devoted to bil- 
liards with circular arcs and, in a smaller extent, to other types of curved boundaries. 
The systematic mathematical description of chaotic billiards and an extended list of refer- 
ences can be found in the book by Chernov and Markarianlsl. Rigorous investigations were 
concentrating on methods for producing fully chaotic billiards and on speciflc properties 
(Bernoulli, K-property, mixing and hyperbolicity) expressing differences between chaotic 



systems 
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iard dynamics have been extensively examined 
16, 17, Q, Q. In recent years, properties of 



classical billiards and their quantum-mechanical counterparts were used to explain and im- 
prove performances of devices in microelectronics and nanotechnolo gy, especially of optical 
microresonators in dielectrical and polymer lasers {20. 21, 22, 23, 24, 25 1. 

We are stressing the fact that notable regions of full chaos have been discovered in bil- 
liards with elliptical arcs and piecewise flat boundaries, indicating that such billiards deserve 



further attention 
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29| . In our previous work we analyzed several types of billiards 



with noncircular arc s (p arabolic, hyperbolic, elliptical and generalized power- law), exhibit- 



ing mixed dynamics 
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321]. Next we investigated, in the full parameter space[33l], the 



elliptical stadium billiards (ESB), flrst introduced by Donnay[7|. Here we extend the same 
type of analysis to the truncated elliptical billiards (TEB), which although similar in appear- 
ance, have different dynamical properties. The truncated elliptical billiard (TEB) is deflned 
by a two-parameter planar domain constructed by truncating an ellipse on opposite sides 
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(Fig. 1). A symmetrical stadium-like shape thus obtained consists of a rectangle with two 
elliptical arcs added at its opposite ends. The corresponding billiard has been introduced 
by Del Magno 29|| who, investigating a restricted part of the parameter space and applying 
the mathematical method of invariant cones, determined the region of hyperbolic behavior 
and presented an estimate of the region where such billiard could be ergodic. 

In the present paper we investigate numerically and analitically the truncated elliptical 
billiard (TEB) in the full parameter space, by using two shape parameters 6 and 7. This 
description of the billiard geometry and dynamics is consistent with our previous analysis of 
the elliptical stadium billiard (ESB) l33l . l34| , which is a two-parameter generalization of 



the Bunimovich stadium billiard ^ and is a special case of the mushroom bi 
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28|, 



liardj35,y,3- 
3^ that this 
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It has been confirmed by analysis and numerical computation j?!. I26l 
billiard is fully chaotic (ergodic) for a sizeable but strictly limited region in the parameter 
space, defined by the stable two-bounce horizontal periodic orbit on one and the pantografic 
orbits on the other side. Our investigations of the ESB and TEB billiards confirm the 
suggestion by Del Magno [29| that in spite of apparently similar stadium-hke shapes, these 
two billiards have essentially different dynamical properties. In the present paper we describe 
our analytical and numerical investigation of the truncated elliptical billiard and compare 
the obtained results with those for the elliptical stadium billiard. 

In Section II we define the TEB billiard boundary and describe its geometrical properties. 
In Section III the existence and stability of selected orbits are discussed and illustrated 
by Poincare plots and orbit diagrams. In Section IV the Poincare sections are used to 
estimate, by means of the box-counting numerical method, the degree of chaoticity for a 
given boundary shape. The results are shown in the parameter-space diagram and compared 
with the same type of diagram for the elliptical stadium billiard. In Section V we briefly 
discuss the possible generalization of the truncated elliptical billiard providing a transition 
between two types of the stadium-like elliptical billiards. Finally, in Section VI we summarize 
the obtained results and propose further investigations. 
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FIG. 1: Three types of the TEB billiard shape: (a) circular shape, with 5 = y^l — 7^: S = 0.6, 
7 = 0.8; (b) flattened shape with elliptical arcs and S > ^J\ — ^'^: 6 = 0.6, 7 = 2.2 and (c) 
elongated shape with S < \/l — 7^: 5 = 0. 6, 7 = 0. 2 



II. GEOMETRICAL PROPERTIES OF THE TRUNCATED ELLIPTICAL BIL- 
LIARD 

In our parametrization the truncated elliptical billiard (TEB) is defined in the x — y plane 
by means of the two parameters S and 7, satisfying conditions < 5 < 1 and < 7 < 00. 
The biUiard boundary is described as 



y{x) 



±7 



,if < |x| < 6 



(1) 



[±7^iEf^ ,iiS<\x\<l 
The horizontal diameter is normalized to 2, so that the horizontal semiaxis of the ellipse is 
1. The vertical semiaxis of the ellipse is 7/vT^^^, and the possible height 27 of the biUiard 
extends from to 00. The horizontal length of the central rectangle is 2S. 

In special cases, for 5 = the shape is a full ellipse, for 5 = 1 it is rectangular, for 
5 = 7 = 1 it is a square and for 5 = and 7 = 1 a full circle. For S = 1/1 — 7^ 
one obtains a set of truncated circle billiards, which separates two distinct billiard classes, 
one for 5 < a/I — 7^ with elongated elliptical arcs and the other with 5 > 1/1 — 7^ and 
fiattened eUiptical arcs. Fig. 1 shows three typical shapes of the truncated elliptical biUiard 
with circular, fiattened and elongated elliptical arcs. 

The coordinates of the focal points are 



for 5 < — 7^, and 



0,± 



(2) 



72 + 52 _ 1 
1-52 



(3) 



for 5 > -y/l — 72. They contain the important term r = 7^ + 5^ — 1 which is negative for 
6 < -y/l — 72, positive for 5 > y^l — 72, and zero for 6 — -y/l — 72 (circular arcs) . This limit 
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FIG. 2: Diagram of the two-dimensional parameter space (7, 5). Lines denote the hmits of existence 
and stabihty for certain orbits, as explained in the text. 



is shown as the thick circular line in Fig. 2 presenting the structure of the 7 — 5 parameter 
space. 

For 5 < |x| < 1 the curvature radius is 

[(i-^^-7^)(i-x^)+7r^ 

7(1 - 5^) 

For < |x| < 6 the boundary is flat and the curvature radius is i? = 00, but for |x| = 5 has 
a discontinuity and drops to 

= ^(1^5^^ 

At the endpoints of the horizontal axis of the ellipse (|a;| = 1) the curvature radius is 

7^ 

^1 = (6) 
which reduces to -Ri = 1 for circular arcs. For full ellipses with 5 = the curvature radius 
at a; = is /2o = 1/7- 

As explained in (33|, the symbols 9, and 0', respectively, denote the angles which the 
normal, the incoming path and the outcoming path make with the x-axis. The angle between 
the incoming (or outcoming) path and the normal to the boundary is /? = (0' — <f))/2. The 
angle between the tangent to the boundary at the point T(x, y) of impact and the incoming 
(or outcoming) path, needed in the computation of the orbit stability, is a = [it/ 2) — (3. 

The angles 6, and 0' are connected by the relation 

2 tan 6 tan + tan 0' 



1 — tan^ 6 1 — tan tan 0' 
The expression (7) is the basis for finding the existence criteria for particular periodic 



orbits[33|]. In our further description we refer to the impact points T{x,y) in the first 
quadrant, with no loss of generality for the obtained results. In the Poincare sections the 
points P(a;, Vx) are obtained by plotting the slope of the velocity direction cos versus 



the x-coordinate of the intersection point with the explained in |30|, l3l|, 122, l33| . 

The Poincare diagrams obtained in this way are area preserving. 
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FIG. 3: Poincare plots for (5 = 0.19 and various 7. 



As described in [13|, l33| , the stability of a periodic orbit is assured if the absolute value 
of the trace of the stability matrix M is smaller than 2, thus if 



- 2 < TrM < 2 



(8) 



Such orbits are elliptic, and those with |TrM| = 2 are neutral (parabolic). The stability 
matrix of the closed orbit of period N can be written as M = M12M23..M7V1, where the 2x2 
matrix Mik for two subsequent impact points Tj and T^, connected by a rectilinear chord 
of the length pjfc, is 



ik 



sm Qj 
sin Qfe 



Pik 



+ 



Pik 



Ri sin ctfc 



Pik 



\ 



sin ai sm 



\ RiRk 



+ 



Ri 



+ 



smofc 



"fc _|_ Pik 
oti Rk sin«i / 



Rk 



(9) 



III. CLASSICAL DYNAMICS OF THE TRUNCATED ELLIPTICAL BILLIARD 



A. Billiards vi^ith b < \J\ 



r 



This subfamily of truncated elliptical billiards has elongated elliptical arcs. In Fig. 3(a-d) 
we show Poincare sections for 5 = 0.19 and different values of 7 < 0.982. Similar results 
for 6 = 0.60 and 7 < 0.80 are shown in Fig. 4(a-d). These pictures reveal a highly chaotic 
behavior. There are no elliptic islands, however, flights of points typical for neutral orbits can 
be discerned. This is remarkably different from the corresponding results for the elliptical 
stadium billiards 33|, where in the same parameter region there were many fixed points and 



elliptic islands due to stable pantographic and other orbits. 



1. The bow-tie orbit 



We investigate the existence and stability of the bow-tie orbit (the lowest pantographic 
orbit), shown in Fig. 5(a). This orbit exists if the coordinates x and y of the impact point 
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FIG. 4: Poincare plots for 6 = 0.60 and various 7. 
and the derivative y' of the boundary at this point satisfy the equation (7), which now reads 



2yy' + x{l - y'^) = (10) 
giving as solution the coordinates of the point of impact 



and 



,1-52-272 ^ ^ 

^= W 1_^2_^2 (11) 



(12) 



The condition 5 < x < 1 that this point should lie on the eUiptical part of the boundary 
leads to the requirement 



This limit is shown in Fig. 2 and is denoted with the letter a. If we denote the points 
with positive x by 1 and the points on the negative side by -1, the deviation matrix can be 
calculated as 

M = {MiiMi_^f (14) 
The corresponding angle a needed in the matrix (9) is given by 



The chords are 



2(1-/- 72) 



272 

p = pi 1 = 2|/ = ^ (16) 

' y(l-52)(l-(52-72) ^ ' 



and 



FIG. 5: Typical lowest order periodic orbits in TEB billiards: (a) bow-tie orbit (neutral); (b) 
rectangular orbit (neutral); (c) horizontal two-bounce orbit (stable); (d) tilted two-bounce orbit 
(neutral); (e) diamond orbit (stable); (f) multidiamond orbit with n = 2 (stable); (g) hour-glass 
orbit (neutral); (h) 8-shaped orbit (stable). 



(J ^ = 27^^ = 2^ (17) 
The curvature radius at the impact point is obtained by substituting (11) into (4) and reads 

(18) 

If we define 

$ = ^ ^_(^ + ^) (19) 

i? sin a i? sin a. R sin a sin a 

the trace of the deviation matrix is 



TrM = 2 [2(2$ + if - l] (20) 

The left-hand side of the stabihty condition (8) is valid automatically, but the right-hand 
side is fulfilled only if 

- 1 < $ < (21) 

By substituting (15), (16), (17) and (18) into (19), one obtains $ = — 1 for all allowed cases. 
The conclusion is that the bow-tie orbit is neutral for all parameter values satisfying the 
existence condition. 



2. The rectangular orbit 

Further we investigate properties of the rectangular orbit shown in Fig. 5(b). According 
to (7), this orbit exists if the derivative on the boundary is y' = — 1. Corresponding solutions 
for the impact point are 



8 



and 



(22) 



^ (23) 

v/(l-52)(l-52+^2) ^ 



The condition 5 < x < 1 leads to the requirement 



6 < 



7 < 



(24) 



2 ' ' ~ 5 

This hmit is shown in Fig. 2 denoted by letter f. Stability is calculated with equation (9) and 
the matrix (14), where the angle a is given by sin a = 1/^/2. The chords are p = pi^i = 2x 
and p' = pi-i = 2y and the curvature radius is 



R 



1-6' 



3/2 



Again, the trace is given by (20), and for this case one obtains $ = 0. The conclusion is 
that also this orbit is neutral for all shapes, both flattened and elongated, allowed by (24). 



29| . Their boundary shapes 



The elongated truncated elliptical billiards were discussed in 
were described by means of two parameters h and a, related to our parameters 6 and 7 as 
follows: 



h = VT^; a= ^^ (26) 

7 



In 



29( 1 the billiards with a > 1 and h < 1 have been analysed and the hyperbolic behavior 
has been identified in the region h < min(l/a, I/V2). In our parameters, this corresponds to 
the quasi-triangular region in the parameter space, denoted by A and B in Fig. 2, delimited 
by curves 6 = ^/T^^ (denoted in Fig. 2 by letter e) and 6 = a/1 — 7^ and by the straight 
line 6 = I/V2. 



The region h < l/vT+~a^ is rigorously proved to be ergodic 29 1. Written with our 
parameters, it obeys the conditions 



72; (27) 



< ^ < ^/i^ (28) 

The corresponding part of the parameter space in Fig. 2 is the one denoted by A. The 
comparison with our results shows that the hmit (27) or (28) is identical to the limit (24) 
in the parameter space, where the parabolic rectangular orbits emerge. 



B. Billiards with 5 > \J\-i^ 

This part of the parameter space, with flattened elliptical arcs, had not been investigated 
previously. In Fig. 3(e-h) we show Poincare sections for b = 0.19 and different values of 
7 > 0.98. In this parameter region dynamics is following the KAM scenario. Similar behavior 
is noticed for values 6 = 0.60 and 7 > 0.89 (Fig. 4(e-h)). Elliptic islands corresponding to 
the horizontal two-bounce and some other orbits are visible, similarly to the corresponding 
results for the elliptical stadium billiards [33]. We investigate the existence and stability 
criteria for these orbits. 



1. Horizontal diametral two-bounce orbits 



The horizontal two-bounce orbit (Fig. 5(c)) obviously exists for all combinations of 6 and 
7, but according to [isl, 33 1 the stability condition 



_P_ 
2R 



< 1 



(29) 



takes the form 



- 7^; 7 > Vl - 5^ (30) 

so that bifurcations giving birth to stable diametral orbits appear at the values S = a/1 — 7^, 
corresponding to circular arcs. In the Poincare diagrams this orbit and the surrounding 
quasiperiodic orbits are visible as two large bands near = ±1. 
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2. Tilted diametral two-bounce orbits 

According to (7) a tilted two-bounce orbit (Fig. 5(d)) exists at the point T{x,y) on the 
biUiard boundary with derivative y' if 

yy' + x = (31) 
This is reahzed for any 6 < x < 1 provided that 

7^ + 5^ = 1 (32) 

thus only for the truncated circle. Since in this case the chord in (29) is p = 2 and the radius 
is it! = 1, these orbits are neutral. 

3. Diamond orbit 

The diamond orbit of period four, shown in Fig. 5(e), exists for any parameter choice. 
It has two bouncing points at the ends of the horizontal semiaxis, and the other two on the 
fiat parts on the boundary. To assess its stability, one should calculate the stability matrix 
M — (MoiMio)^. The angles contained in the matrix are given as 



7 . 1 

sm CKo = — ; sm ai — - 

P P 



where 



(33) 



p = + 7' 

The curvature radius at x = 1 is given by (6). This leads to the trace 



TrM = 2 



(34) 



(35) 



and to the condition < i? or 1 + 7^ < 7^/(1 — 5^), thus the stable diamond orbit appears 
when 
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This limit is shown in Fig. 2 as the hne denoted by letter c. 



4. Multidiamond orbits 

The multidiamond orbit of order n is the orbit of period 2 + 2n, which has two bouncing 
points at the ends of the horizontal axis and 2n bouncing points on the flat parts of the 
boundary (Fig. 4 (f)). Such orbit exists if 



6 > 1 



1 



n 



(37) 



As explained for a similar case in 331], the chord p in (35) should be replaced by 



L = npr, 



(38) 



where, for the truncated elliptical billiard. 



(39) 



The trace of the stability matrix is then 



TrM = 2 



(40) 



with R given by (6). The resulting condition for the stability of the multidiamond orbit is 



5>Jl- ^ ^\ 7>a/ , (41) 

The limiting curves (41) are plotted in Fig. 2. The line with n = 2 is denoted by letter d, 
and above it there are several lines for n > 2. For 7 ^ 00 the minimal values of 6 above 
which the multidiamond orbits appear are 
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FIG. 6: Poincare plots for a set of different values (5 = 7, showing the appearance of successive 
multidiamond orbits of higher order. Two symmetrical triangle-shaped islands for (5 = 7 = 0.90 
originate from the "8-shaped" orbit. 



7-+00 V ffi 



\iYa5 = \ll-\ (42) 



The emergence of multidiamond orbits can be followed by observing the Poincare sections 
for a set of shapes with 5 = 7 (Fig. 6). The values of this parameter for which an orbit 

of new n appears are given as intersections of the straight line 5 = 7 with curves (41), and 
obey the equation 

n'^S^ - (n^ - 2)52 -1 = (43) 
For the diamond orbit (n = 1) this equation reads 

54 + 52 - 1 = (44) 

and the orbit appears for 



5 = 7= y(\/5- l)/2 = 0.78615. 
For the same type of boundary the stable two-bounce orbit appeared at 

5 = 7 = I/V2 = 0.707107. 



5. The hour-glass orbit 

The hour-glass orbit (Fig. 5(g)) looks hke the bow-tie orbit rotated by 7r/2. It exists if 
the coordinates x and y of the impact point and the derivative y' of the boundary at this 
point satisfy the equation 

2xy' + yy'^ - y = (45) 
giving as solution the coordinates of the impact point 
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and 



(46) 



The condition 6 < x < 1 that this point should he on the eUiptical part of the boundary 
leads to the requirement 



V2(l - 5^) < 7 < (48) 

These limits define the region shown in Fig. 2 denoted by letter b. 

If we denote the points with positive y by 1 and the points on the negative side by -1, 
the deviation matrix can be calculated from (14). The angle a needed in the calculation is 
given as 



^^"^=^/ 2(^^77^-l) 



The curvature radius at this point is given as 



The chords are 



P-Pi,i = 2. = 2y^^l-|^ (51) 
and 

p' ^ = 2v/^^ = ^^^^ 
If we define $ as in (19), the trace of the deviation matrix is again given by (20) and the 
orbit is stable if -1 < $ < 0. 
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FIG. 7: Diagram showing the degree of chaoticity gdass of the truncated ehiptical bilhards (TEB), 
in dependence on the shape parameters, with < 7 < 2.2 and < 6 < 1. Black points denote 
shapes with gdass = 1-00, red with 0.90 < Qdass < 1-00, green with 0.80 < Qdass < 0.90, yehow with 
0.70 < gciass < 0.80, blue with 0.60 < gciass < 0.70 and grey with 0.00 < gdass < 0.60. 

When we substitute the calculated values of R, p, p' and sin a into (19), we obtain $ = — 1 
for all allowed shapes and conclude that the hour-glass orbit is neutral. This means that in 
the truncated elliptical billiards (TEB) there is no stable hour-glass orbit, at variance with 
the elliptical stadium billiard (ESB), where such an orbit having interesting properties was 
stable in a large fraction of the parameter space 33|] . Besides the diamond and multidiamond 
orbits, in Fig. 6 one discerns the presence of another, "8-shaped", stable orbit, shown in 
Fig. 5(h). 

IV. THE BOX-COUNTING NUMERICAL ANALYSIS OF THE DEGREE OF 
CHAOTICITY IN THE FULL PARAMETER SPACE 

In this section we return to the question of limits within which the truncated elliptical 
billiard is fully chaotic. Here we test these limits numerically, with the help of the box- 



counting method[l4, l33|, l38|. We calculate the Poincare sections for a chosen pair of shape 
parameters, starting with n\ randomly chosen sets of initial conditions and iterating each 
orbit for ^2 intersections with the x-axis, thus obtaining n\ x n2 points in the Poincare 
diagram. Then we divide the first quadrant of the phase plane into a grid of n x squares 
(boxes), count the number of boxes which have points in them and calculate the ratio of 
this number to the total number of boxes. The obtained ratio is denoted by gdass- In this 
way also certain points belonging to invariant curves within the regular islands are included. 
But since our main aim is to examine the onset of full chaos, this method gives satisfactory 
results, providing that the appropriate values of rii, and n are used. Detailed testing has 
shown that reliable results are obtained for values n\ = 100, 77.2 = 5000 and n = 100 used in 
our present calculation [J 



In Fig. 7 we plot in the 5 — 7 plane the points representing the pairs of shape parameters. 
Points are plotted in different colors, depending on the corresponding value of gdass- The 
full chaos, corresponding to gciass = 1.00, is depicted by black points. Colored points denote 
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FIG. 8: Diagram showing the degree of chaoticity ^class of the ehiptical stadium bilhards (ESB), 
in dependence on the shape parameters, with < 7 < 2.2 and < 6 < 1. Black points denote 
shapes with gciass = 1-00, red with 0.90 < gciass < 1-00, green with 0.80 < gciass < 0.90, yellow with 
0.70 < Qciass < 0.80, blue with 0.60 < ^ciass < 0.70 and grey with 0.00 < ^ciass < 0.60. 

FIG. 9: Gonstruction and parameters of the generalized elliptical truncated stadium-like billiards 
(GTESB) 

shapes within intervals between and 0.99. This diagram confirms that for the truncated 
elliptical billiards (TEB), in the region below the onset of the stable two-bounce horizontal 
orbit, dynamics is practically completely chaotic. This is in strong contrast with the behavior 
of the elliptical stadium billiard for which the similar diagram is shown in Fig. 8. For the ESB 
billiards the region of chaos was strictly bounded also from the lower side and determined 
by emergence of stable pantographic orbits. 

To examine the possible mechanism for this difference, we assume that the TEB and the 
ESB billiards are two extreme cases and search for a possible transition between them. 

V. GENERALIZED TRUNCATED STADIUM-LIKE ELLIPTICAL BILLIARDS 

In this section we propose a new large class of stadium- like billiards which we call gener- 
alized truncated elliptical stadium-like billiards (GTESB). Such a billiard depends on three 
shape parameters 5, 7 and k. The allowed values of the shape parameters are 



Q<S <1; 0<7<oo; -1<«;<1 (53) 

For the limiting values of k we obtain the two billiard families considered before: ior k — — 1 
GTESB reduces to the elliptical stadium billiard (ESB), and for k = 1 GTESB becomes the 
truncated elliptical billiard (TEB). 

The new GTESB billiard boundary is obtained by adding elliptical arcs symmetrically at 
the two opposite ends of a rectangle with sides 2S and 2j. Elliptical arcs are cut out from 
the two identical but generally detached ellipses by two horizontal straight lines at y = ±7 



16 



FIG. 10: Poincare plots for the generalized truncated elliptical stadium-like billiards (GTESB), for 
6 = 0.2 and 7 = 0.6 with different values of k: (a) k = —1.0; (b) k = —0.6; (c) k = —0.3; (d) 
K = 0;(e) K = 0.3; (f) k = 0.6; (g) k = 0.90; (h) k = 1.0. 

(Fig. 9). The two ellipses have centers at the points 

X± = ±S (^^) (54) 

The distance between the two centers is D — 5(1 — k). The horizontal and vertical semiaxis 
are given, respectively, as 

A.^l-S (^i^^ (55) 

and 



7 



and the equation of the two ellipses reads 



S 



1 



(56) 



The horizontal diameter of the billiard is 2. For k — —1 and ^ — 1 — S the GTESB becomes 
the Bunimovich stadium billiard. 

In Fig. 10 the Poincare sections are shown for 6 = 0.2 and 7 = 0.6 with k assuming 
different values between -1 and 1. The four islands typical for the bow-tie orbit arc present 
for all K except for k — 1 (TEB), where this orbit becomes neutral and the island reduces 
to a caracteristical flight of points. The limits separating chaotic from mixed behavior are 
determined by the onset of the stable horizontal 2-bounce orbit and are given by (29). Since 
the curvature radius at |x| = 1 is 

_ 7^[2-^(l-/.)] 

- 2{1-S){1 + 5k) ^^^^ 
the upper limit of chaos is determined by the condition 
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FIG. 11: Dependence of gciass on 7 for a set of shapes with (5 = 7 and the choice of k as in Fig.lO. 



In Fig. 11 the chaotic fraction Qdass is shown for the special case 5 = 7 for different k, in 
dependence on 7. It is noticed that in the case k — —1 (ESB) there is a narrow, strictly 
limited region of full chaos, outside of which the values of gciass are low. For k, = 1 (TEB) 
the fully chaotic region is much larger and extends practically over all values S and 7 below 
the chaotic limit. Between these two limits, the regions of full chaos (gciass — 1) are shorter 
and limited, but there are many shapes with chaotic parameter close to 1 (between 0.90 and 
0.99). This corresponds to a selection of narrow islands in the Poincare plots, as seen in Fig. 



In conclusion, our investigation of the elliptical stadium-like billiards has revealed a rich 
variety of integrable, mixed and chaotic behavior, which is connected with the character of 
the two elliptical arcs and with their mutual position. This strong dependence on parameters 
5 and 7 is confirmed for the truncated elliptical billiards, but is even more enhanced when 
a third shape parameter k is added. Analysis shows, however, that among all considered 
shapes the truncated elliptical billiard (TEB), created by cutting a single ellipse with two 
parallel straight lines, has exceptional properties, notably that it is chaotic practically in the 
whole region of elongated elliptical arcs. Notable is the presence of many neutral orbits in 
this region, consistent with the fact that these orbits actually can be identified as orbits in 
an ellipse. For the flattened arcs, the stable islands due to the two-bounce horizontal orbit 
and to the diamond and multidiamond orbits occupy an important part of the phase plane. 

Our investigations can be useful for the experimental application of billiards in the laser 
technology, where properties and directional intensities of the optical microresonators depend 
strongly on the boundary shape. They can also be applied in designing the semiconducting 
optical devices and in the technology of microwave and acoustic resonant cavities. With this 
purpose in mind, we propose further analysis of the stadium-like billiards with elliptical arcs 
and an extension of the present investigation to different types of open billiards. 




(59) 
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